Critical Point of an Interacting Two-Dimensional Atomic Bose Gas 
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We have measured the critical atom number in an array of harmonically trapped two-dimensional 
(2D) Bose gases of rubidium atoms at different temperatures. We found this number to be about 
five times higher than predicted by the semi-classical theory of Bose- Einstein condensation (BEC) in 
the ideal gas. This demonstrates that the conventional BEC picture is inapplicable in an interacting 
2D atomic gas, in sharp contrast to the three-dimensional case. A simple heuristic model based on 
the Berezinskii-Kosterlitz-Thouless theory of 2D superfluidity and the local density approximation 
accounts well for our experimental results. 

PACS numbers: 03.75.Lm, 32.80.Pj, 67.40.-w 
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Bose-Einstein condensation (BEC) at a finite temper- 
ature is not possible in a homogeneous two-dimensional 
(2D) system, but an interacting Bose fluid can neverthe- 
less become superfiuid at a finite critical temperature [H . 
This unconventional phase transition is described by the 
Berezinskii-Kosterlitz-Thouless (BKT) theory and 
does not involve any spontaneous symmetry breaking and 
emergence of a uniform order parameter. It is instead as- 
sociated with a topological order embodied in the pairing 
of vortices with opposite circulations; true long-range or- 
der is destroyed by long wavelength phase fluctuations 
even in the superfiuid state [1, [B| . 

Recent advances in producing harmonically trapped, 
weakly interacting (quasi-) 2D atomic gases [g, 0,11, @, 
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[12j, LL3J, LLJ] have opened the possibility for de- 
tailed studies of BKT physics in a controllable envi- 
ronment. There has been some theoretical debate on 
the nature of the superfiuid transition in these sys- 
tems 15, M E3, d m 

because the harmonic con- 
finement modifies the density of states compared to the 
homogenous case. This allows for "conventional" finite 
tem per ature Bose-Einstein condensation in the ideal 2D 
gas [20| . Early experiments have been equally consistent 
with the BEC and the BKT picture of the phase transi- 
tion. For example, the density profiles at very low tem- 
peratures @ are expected to be the same in both cases. 
However, recent studies of matter wave interference of in- 
dependent 2D atomic clouds close to the transition have 
revealed both thermally activated vortices 12, [l3[ and 
quasi-long-range coherenceproperties [l3| in agreement 
with the BKT theory dHHl]. 

In this Letter, we study the critical atom number in an 
array of 2D gases of rubidium atoms, and observe stark 
disagreement with the predictions of the ideal gas BEC 
theory. We detect the critical point by measuring (i) the 
onset of bimodality in the atomic density distribution 
and (ii) the onset of interference between independent 
2D clouds. These two measurements agree with each 
other, and for the investigated range of temperatures 
T « 50-110nK give critical atom numbers N c which 
are ~ 5 times higher than the ideal gas prediction for 
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Figure 1: Experimental setup. Left: A one-dimensional 
optical lattice is used to split a magnetically trapped 3D 
BEC (transparent ellipsoid) into a small array of 2D clouds. 
Right: Contour lines of the total (magnetic and light) poten- 
tial V(x, z) in the y = plane for the lattice phase such that 
the two central planes are symmetric with respect to the trap 
center. 

conventional Bose-Einstein condensation in our trap [20l ] . 
For comparison, in three-dimensional (3D) atomic gases, 
where conventional BEC occurs, the increase of the criti- 
cal atom number due to repulsive interactions is typically 
on the order of ten percent [H, [24J. A simple heuris- 
tic model based on the BKT theory of 2D superfluidity 
and the local density approximation gives good agree- 
ment with our measurements. 

In [l3[ we studied quasi-long-range coherence of a 
trapped 2D gas,which is directly related to the super- 
fluid density p s [22j . In that case, signatures of the BKT 
transition emerge only once a significant part of the cloud 
becomes superfiuid. Since the atomic density in the trap 
is not uniform, this happens slightly below the true criti- 
cal temperature for the onset of superfluidity in the trap 
center, and the observed transition is rounded off. The 
present study concentrates on the exact critical point and 
relates to the total density at criticality p c , which has 



25, 2 



been of long standing theoretical interest 

Our experimental procedure for the preparation of cold 
2D Bose gases has been described in [131] . We start with 
a 87 Rb 3D condensate in a cylindrically symmetric mag- 
netic trap with trapping frequencies lo x = 2n x 10.6 Hz 
and L0y — lu z = 2n x 125 Hz. To split the sample into 
2D clouds we add a blue detuned one-dimensional opti- 
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Figure 2: Phase transition in a rubidium 2D gas. 2D clouds 
confined parallel to the xy plane are released from an optical 
lattice and the density distribution is recorded by absorption 
imaging along y after t — 22 ms of time of flight. The mea- 
sured line densities p(x) (•) for an atom number just below 
(left) and just above (right) the critical number are displayed 
together with bimodal fits (solid lines) . The dashed line in the 
left panel shows the expected distribution of the 2D ideal gas 
at the threshold of conventional BEC in our potential at the 
same temperature (T = 92 nK). The dotted line in the right 
panel indicates the Gaussian part of the bimodal distribution. 

cal lattice with a period of d = 3 /im along the vertical 
direction z (see Fig. [TJ). The lattice is formed by two 
laser beams with a 532 nm wavelength and focussed to 
waists of about 120 /im, which propagate in the yz plane 
and intersect at a small angle. The depth of the lattice 
potential around x = is h x 35 kHz, corresponding to 
a vertical confinement of tu z — 2tt x 3.0 kHz. The tun- 
neling rate between adjacent sites at the center of the 
trap (x = 0) is negligible on the time scale of the exper- 
iment. The finite waists of the lattice beams result in a 
slow variation of u z along x, and the variation of the zero 
point energy Hu> z (x)/2 modifies ui x to 2ir x 9.4 Hz at the 
trap center. 

Fig. Q] shows contour lines for the full trapping po- 
tential. The number of significantly populated lattice 
planes is ~ 2 — 4 in the investigated temperature range 
(50-110 nK). The vast majority of atoms is trapped in 
the central x region where the 2D criterion kT < fiio z (x) 
is fulfilled. However, the exchange of particles between 
lattice sites is still possible via the far wings of the en- 
ergy distribution (at energies above 460 nK). This en- 
sures thermal equilibrium between the planes [28j on the 
time scale of ~ 100 collision times |29j . which in our case 
corresponds to a fraction of a second. The 2D interac- 
tion strength is g — (h 2 /m)g, where the dimensionless 
coupling constant g = a s \J% r Kmuj z lh — 0.13, a s = 5.2 nm 



is the scattering length, and m the atomic mass [15l. |27|. 
The interaction energy Ei nt ~ gpo, where po is the peak 
density, also satisfies the 2D criterion E- lrA < hu) z . 

We measure the critical atom number N c by varying 
the total atom number N at a fixed temperature. We 
start with a highly degenerate sample and keep it trapped 
for a time r varying between 1 and 10 s. During this 
time we maintain a constant temperature by applying a 
constant radio frequency field in the range of 10-25 kHz 



above the frequency corresponding to the bottom of the 
trap. As the hold time r increases, N gradually reduces 
and drops below N c due to inelastic losses. 

The atomic density profiles are recorded in the xz 
plane by resonant absorption imaging along y after t = 
22 ms of time of flight expansion. Along z the profiles are 
Gaussian, closely corresponding to the zero point kinetic 
energy fku z (x — 0)/4. Along x, for all N < N c a Gaus- 
sian distribution fits the data well [3(|. For N > N c , the 
profiles exhibit a clearly bimodal shape (Fig. [2]). The bi- 
modal distributions are fitted well by the sum of a Gaus- 
sian, corresponding to the "normal component" , and a 
parabolic Thomas-Fermi (TF) profile expected from su- 
perfluid hydrodynamics [231 ] . 

From the bimodal fits we extract the total atom num- 
ber N and the number of atoms within the TF part of 
the distribution N . The absolute detection efficiency of 
our imaging system was calibrated by measuring critical 
atom numbers for 3D BECs, taking into account interac- 
tion effects 23|, 24 1 . For a given energy of the evaporation 



surface E cvlip the width of the Gaussian part of the dis- 
tribution is nearly independent of N (see inset of Fig. [3]). 
For a quasi-non degenerate gas (N ~ N c /2) this width 
is given by the temperature. We thus use this estimate 
for T also in the degenerate regime (N > N c ), although 
one could have expected in this regime a deviation from 
the Gaussian law for the normal fraction. For example 
the ideal gas theory predicts a distribution at the BEC 
point that is much more peaked at the center of the cloud 
(Fig. [2]). The temperatures inferred in this way scale as 
Ecv&p — V^T, with rj ss 10 compatible with the usual 3D 
values for evaporation times equal to a few hundred colli- 
sion times [29( . We estimate the systematic uncertainties 
of our atom number and temperature calibrations to be 
20% and 10%, respectively. 

Fig. [3] illustrates the threshold behavior of Nq j3lj, and 
Fig. [4] shows the critical numbers N c measured at four 
different tem per atures. In a single 2D ideal gas, BEC is 
expected for 
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where w is the geometric mean of the two trapping fre- 
quencies in the plane. For comparison with our experi- 
mental results we have numerically integrated the Bose- 
Einstein distribution for our confining potential sketched 
in Fig. [TJ The result depends on a ten percent level on 
the exact position of the lattice planes relative to the 
minimum of the magnetic trap potential. Since we do 
not fully control this position, we average over the possi- 
ble configurations. We obtain the result N™^ u — pN Ct id, 
where the effective number of planes p smoothly grows 
from w 2.2 at 50 nK to w 4.2 at 110 nK. The resulting 
N™^ U (T) is shown in Fig. [4] as a solid line. Our mea- 
surements clearly show systematically higher N c than ex- 
pected for ideal gas condensation. An empirical function 
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Figure 3: Measurement of the critical point. The number of 
atoms in the Thomas-Fermi part of the bimodal distribution 
iVo (o) is plotted as a function of the total atom number N. 
The solid line shows the linear fit we use to determine N c , and 
the dashed line is its extrapolation. For comparison, the inter- 
ference amplitude |/3(0, fco)| (•) is also displayed as a function 
of N. It shows the same threshold N c within our experimen- 
tal precision. The inset shows that the temperature deduced 
from the Gaussian part of the fit is to a good approximation 
constant for all data points. Horizontal and vertical dashed 
lines indicate the average temperature and the critical atom 
number, respectively. The solid line marks the region used to 
determine the average temperature T = 92(6) nK close to the 
transition. Each data point is based on 5-10 images, all error 
bars represent standard deviations. 
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Figure 4: Critical point in an interacting 2D gas. The crit- 
ical atom number iV c (•) is measured at four different tem- 
peratures T. Displayed error bars are statistical. The solid 
line shows the ideal 2D gas BEC prediction N™$ u . The 
dashed line is the best empirical fit to the data, which gives 
N c =aN^ u with a = 5.3(5). 

iV c = aN™$ ti (T), with the scaling factor a as the only 
free parameter, fits the data well and gives a — 5.3(5), 
where the quoted error is statistical. 

We also study the coherent fraction of the 2D gas and 
compare its behavior with the bimodal density profiles. 
We investigate the interference patterns that form af- 
ter releasing the independent planar gases from the trap 
(Fig. [13j. Fourier transforming the density profile 
p{x, z) — > T\p{x, z)\ = p(x, kg) allows us to quantify the 
size of the coherent, i.e. interfering part of the gas as 
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Figure 5: Coherence and density profile below the transition. 
Interference of 2D clouds is used to compare the coherent 
part of the cloud with the part following the Thomas-Fermi 
density distribution. Left: Interference patterns in the xz 
plane (see example in inset) are Fourier transformed along 
the expansion axis z and averaged over ten images taken un- 
der identical conditions, to obtain \p(x,k z )\. Right: Within 
experimental precision, fits to the total density profile \p(x, 0)| 
and the interference amplitude profile \p(x, ko)\ give the same 
Thomas- Fermi diameter 2J?tp , indicated by the dashed lines. 
The weak second harmonic peak at k z — 2ko reveals small 
occupation of the outer lattice planes. 

a function of N. The spatial frequency corresponding 
to the fringe period for the interference of neighboring 
planes is ko = md/ht. We find that p(x 1 ko) is well fitted 
by a pure Thomas-Fermi profile. Within our experimen- 
tal accuracy, the radii -Rtf(^o) of these profiles are equal 
to those obtained from a bimodal fit to the density. In 
particular, the onsets of interference and bimodality co- 
incide (circles and diamonds in Fig. [31 respectively) . 

We now turn to the interpretation of our measurements 
in the framework of the BKT theory of 2D superfluidity. 
The theory predicts a universal jump of the superfluid 
density at the transition, from p s — to p s X 2 = 4, where 
A = hj\/2nmkT is the thermal wavelength [2l[ (for ex- 
periments see [l|, HH). However, the total density at the 
critical point p c is not universal because it depends on the 
microscopic interactions. For weak interactions (g < 1), 
p c A 2 = \n{C/g) [H, with C = 380 ± 3 given by high- 
precision Monte Carlo calculations (26[. For our value of 
g = 0.13 (experimentally confirmed by measuring -Rtf 
as a function of No) this gives p c X 2 = 8.0. 

In a harmonic trap, within the local density approxi- 
mation, the transition is expected to occur when the den- 
sity in the center of the cloud reaches the critical value 
p c . We can heuristically relate the critical density and 
the corresponding critical atom number 2V C bkt using the 
experimentally observed Gaussian density profiles. For a 
single plane with a quadratic confining potential V(x, y), 
integrating p(x, y) = p c exp(- V(x, y)/ (kT)) gives: 

, fkT\ 2 9 6 

iV c ,BKT = Pc A 2 f — J = p c A 2 iV c ,id • (2) 

For p c A 2 = 8.0 this gives iV^BKT/^Vcad = 4.9. In a lattice 
configuration this ratio changes only slightly. We set the 
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peak density in the most populated plane to p c , and sum 
the contributions of all planes j using the corresponding 
potentials Vj (here we neglect the small non-harmonic 
effects due to finite laser waists). The total population in 
the lattice is then iV^^T = p'-/V Cj bkt, where the effective 
number of planes p' varies from 2.4 at 50 nK to 3.5 at 
110 nK. We thus obtain N^^ T /N^ ti ~ 4.7, which is 
close to the experimental ratio a = 5.3. 

One could try to reproduce our observations within the 
self-consistent Hartree-Fock (HF) theory [h| (see also 
H, H, 34j|). by replacing V(r) with the effective mean 
field potential V(r) + 2gp(r) and again setting the peak 
density to the BKT threshold p c . For very weak interac- 
tions, log(l/<?) » 1, analytical HF calculation gives criti- 
cal numbers which are only slightly larger than iV Ci id [lij ]. 
This approach could in principle be implemented numer- 
ically for our value of g and our lattice geometry. How- 
ever, it has been suggested [35| that treating interactions 
at the mean field level is insufficient for g ~ lO^ 1 , be- 
cause the interactions are strong enough for the critical 
region to be a significant fraction of the sample. In fu- 
ture experiments with atomic gases g could be varied 
between 1 and 10 -4 using Feshbach resonances, allowing 
for detailed tests of the microscopic BKT theory and the 
possible breakdown of the mean field approximation. 

In conclusion, we have shown that the ideal gas the- 
ory of Bose-Einstein condensation, which is extremely 
successful in 3D, cannot be used to predict the critical 
point in interacting 2D atomic gases, where interactions 
play a profound role even in the normal state. A much 
better prediction of the critical point is provided by the 
BKT theory of 2D superfluidity. We have also shown 
that, despite the absence of true long-range order, the 
low temperature state displays density profiles and local 
coherence largely analogous to 3D BECs. 
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